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Abstract 

We determine the 0{as) radiative corrections to polarized top quark pair production 
in e~^e~ annihilations with a specified gluon energy cut. We write down fully analyt- 
ical results for the unpolarized and polarized 0{as) cross sections e^e~ — > it{G) and 
e'^e~ — > W (G) including their polar orientation dependence relative to the beam 
direction. In the soft-gluon limit we recover the usual factorizing form known from 
the soft-gluon approximation. In the limit when the gluon energy cut takes its max- 
imum value we recover the totally inclusive unpolarized and polarized cross sections 
calculated previously. We provide some numerical results on the cut-off dependence 
of the various polarized and unpolarized cross sections and discuss how the exact 
results numerically differ from the approximate soft-gluon results. 



1 Introduction 



After the discovery of the heavy top quark at the Tevatron in 1995 there has been much 
interest in the use of the proposed high energy hnear e~^e~ colhder as a copious source 
of top quark pairs. When the proposed hnear colhder ILC comes into operation it is 
necessary to have available detailed radiative corrections to the production and the decay 
of top quark pairs. Concerning production there are a number of unpolarized and single 
spin polarized structure functions that describe the e~^e~ production process of massive top 
quark pairs. In the unpolarized case one has the three structure functions Hu (unpolarized 
transverse), (longitudinal), and Hp (forward-backward) which determine the polar 
angle orientation of the top pair relative to the beam axis. Partial results on the full 
0{as) radiative corrections to the unpolarized structure functions Hu, and Hp had 
been written down in Refs. [HE] starting with the early work on the 0{a) QED radiative 
corrections to the vector current (7ve"'"e~) vertex function Complete results on the 
0{as) unpolarized structure functions have been first given in Refs. [H [5]. All of the 
unpolarized 0{as) structure functions were recalculated in the course of computing the 
top quark's 0{as) polarization asymmetries where the unpolarized structure functions 
were needed to normalize the polarization asymmetries [6l [TJ El [HI [ID] • The numerators of 
the polarization asymmetries are expressed in terms of polarized structure functions. In 
the case of the longitudinal polarization of the top one has the three structure functions 
H\j, H'^ and Hp for which the full 0{as) radiative corrections were given in Refs. [7t[8|[T0]. 
In the case of a top quark polarized transverse or normal to the event plane, one has two 
structure functions in each case which are Hj and ifj, and Hf^ and H^, respectively (see 
e.g. Ref. [9]). These were calculated in Refs. P ITT]. 

When doing the full 0{as) radiative corrections one integrates over the full (hard and 
soft) gluon phase space. For some applications it is also interesting to consider radiative 
corrections where one integrates over gluon phase space up to a given gluon energy cut 
Such radiative corrections may be dictated by experimental considerations when soft 
gluons accompanying the top quark pair cannot be resolved by the detector. Alternatively 
one could attempt to measure the cross section for top-antitop-gluon production with a 
given gluon energy cut Ec and compare the energy cut dependence of the cross section 
with the predictions of QCD. Finally, one could define a hard gluon region by introducing 
a lower gluon energy cut and compare experiment with QCD in the hard gluon region. 

In this paper we provide analytical results for the 0{as) radiative corrections to the 
three unpolarized structure functions Hu, H^, and Hp as well as for the seven polarized 
structure functions H[j, Hf^, Hp, hJ'^ , and -ffj'^ for polarized top quarks where we 
integrate over the gluon energy phase space up to a given energy cut Ec. We mention that 
radiative corrections with a gluon energy cut have been treated before in the unpolarized 
case [I21[I3]. 

We emphasize that we are not using the soft-gluon approximation (SGA) in the present 
calculation but integrate over the full 0{as) matrix element tree graph structure. However, 

^Technically, this means that one is dealing with a three-scale problem with the scales q^, nit and Ec- 
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we will compare our results with the soft-gluon approximation. The soft-gluon approxi- 
mation consists of the factorization of the tree graph contribution into the Born term 
contribution and a universal soft-gluon piece which can be easily integrated. An 0{as) 
calculation of some of the structure functions appearing in polarized top pair production 
using variants of the soft-gluon approximation has been done before in Refs. [TOl IT^ . 

One of the further aims of the present investigation is to find out to what extent one can 
pin down a new non-SM (Standard Model) coupling structure in top quark pair production 
in the presence of 0{as) radiative corrections with an exact treatment of gluon emission 
rather than soft-gluon emission. In the latter approximation the tree graph contribution 
is Born term like and thus polarization-type observables would not be affected by the ra- 
diative tree graph corrections but only by the non-Born term structure of the one-loop 
contributions. Deviations from SM predictions for the polarization-type observables could 
result from new non-SM coupling structure or from an exact treatment of radiative correc- 
tions. As an example we will introduce an anomalous CP-odd axial current and compare 
the results of our exact next-to-leading order (NLO) calculation with the contributions of 
the anomalous axial current for some relevant observables and structure functions. 



2 Unpolarized and polarized structure functions 

In order to acquaint the reader with our notation, we use this section to outline the main 
structure of the cross section calculation and to introduce the various unpolarized and 
polarized structure functions that come into play. To start with, we define a polarized 
hadron tensor for the three-body process (7^, Z) qijpi) + ^(^2) + G{pz) according to 

H,M.Pi.P2.s)= m-'^)G\jMml\qq{s)G) (1) 

g,G spins 

where pi, p2 and ps are the four- momenta of the quark, antiquark and gluon, respectively, 
and 5' = pi + + is the four-momentum of the intermediate gauge boson. The spin 
vector of the quark is denoted by s. A similar definition holds for the Born case (7^, Z) 
q{pi) + q{p2)- The hadron tensor defined in Eq. ([T]) depends on the vector (V^: 7^) and 
axial- vector {A: 7^75) composition of the product of currents and jy. It is convenient 
to introduce the four independent hadron tensor components iJ^^, {i = 1,2,3,4) defined 
according to 

tt3 _ '^(TjVA ttAVn tt4 _ ^ ( tt'^A . ttAV\ (r,\ 

For notational convenience we have omitted all arguments in the hadron tensor components 
in Eqs. ([2]). In the following we will use explicit arguments only when they are needed. For 
example, we include the spin vector argument when we define unpolarized and polarized 
structure functions if*^ and -ff^'™ {i = 1, 2, 3,4, m = i, T, N) according to 

H' = HUsn + HU-sn, Hjir = HUsn - HU-s^ (S) 
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Figure 1: Orthonormal spin basis er, gn and for the top quark. Also shown are the 
beam plane (light gray, respectively, yellow) and the event plane (dark gray, respectively, 
green) . 



where is the spin vector corresponding to longitudinal (m = £), transverse (m = T) 
and normal (m = A^) polarization of the top quark. Our choices of the three orthonormal 
spin directions (ey, e^, e^) are given by 

{Pe- X pi) xpi ^ pe- xpi ^ Pi ... 

er = cn = = — — (4) 

\{Pe- X pi) X pi\ X pi\ 

(cf. Fig. [1]). For the hadron tensor components we introduce the compact notation H^J;^^ 
where the round brackets indicate that, in the unpolarized case, the index m and the round 
bracket is omitted. We use this compact notation to display the general features common 
to the unpolarized and polarized parts. 

For the process e^e~ QQiG), the cross section can be written in modular form 
consisting of the hadron tensor, the lepton tensor and the model dependent coupling co- 
efficients gij. The Standard Model (SM) values of the coupling coefficients gij are listed in 
Appendix A. The unpolarized and polarized cross sections read 

^^^'"^ = 1^ E g^,L''''H^^r)dPS (5) 

where dPS is the phase-space factor. The lepton tensor components L^^ {i = 1,2,3,4) 
are defined in the same way as in Eq. ([2]). The process e~^e~ (iq{G) can be described 
either in the beam plane spanned by the electron and positron beam and the outgoing 
quark, or the event plane spanned by the quark, the antiquark, and the gluon. In the 
Born case where no gluon is emitted, both planes coincide by convention. The polar angle 
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between the quark momentum and the electron momentum is denoted by 9 (or by Ote-)-, 
and the azimuthal angle between the two planes is denoted by x- order to determine 
directions, we define different frames with the (x, z) plane lying in the corresponding plane. 
For the beam plane we define a lepton frame with the z axis determined by the momentum 
direction of the electron, and a beam frame with the z axis determined by the momentum 
direction of the quark. For the event plane we define an event frame with the z direction 
determined again by the momentum direction of the quark. The transition from one frame 
to the other is performed by using the two Euler angles 6 and x- 

The natural frame for describing the hadron tensor is the event frame which makes 
no reference to the beam plane. On the other hand, the lepton tensor is most naturally 
described in the lepton frame. In this frame the lepton tensor component L^t^" vanishes 
identically and L^^'^ vanishes for zero lepton masses (which we assume). The remaining 
two components have the simple form 
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(6) 



The contraction of the lepton and hadron tensor has to be done in one particular frame 
for which we choose the event frame. We therefore have to rotate the lepton tensor into 
the event frame. In doing so a variety of angular dependences appear. In fact we can 
decompose the lepton tensors according to 

L^'^" = ^ |i(l + cos' 0)77^ + sin^ ^ iir -2V2 sin ^ cos ^iir}, 
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{cos^iir -2\/2sin^i77} 



(7) 



where Uj and Ua contain an implicit linear dependence on sin x and cos x- The matrices 
77(7, Uj, Up and Ua are called projectors because when contracting the lepton tensor 
with the hadron tensor they project out the relevant coefficients of the hadron tensor that 
give rise to the various angular dependences. The decomposition in Eq. ([7]) describes the 
complete angular dependence of unpolarized and polarized top production in the process 
e~^e~ qqiG). It gives rise to the decomposition of the differential cross section according 
to 



d cos 9 8 
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(9) 



Without beam polarization effects one finds the following pattern. For i = 1 one has 



contributions from a 



U, L, I and for i 



4 one has contributions from a = F, A as 
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written out in Eq. ([7]). More details about the coupling pattern including transverse and 
longitudinal beam polarization effects can be found in [10|. In Eqs. ([9]) we have divided out 
the d cos 6 differential which has already been taken into account in the polar distribution 
(IH]). For the two-particle final states (Born term and loop contribution) one has the phase- 
space factor 

dPS2 = dcosOdx^ — dcos9. (10) 

8(27r)2 ^ 167r ^ ' 



where v = yl — Arn^/q^ is the velocity of the outgoing quark. The transition to the 
rightmost form in Eq. ( fTOj) marked by an arrow expresses the fact that the azimuthal 
integration over x is always implied throughout this paper. As we shall see, the transverse 
and normal spin dependence drop out for the components a = U, L, F in Hl^"^'' but are 
retained for the components a = I,A after the azimuthal integration over Just the 
opposite happens to the spin independent and longitudinal spin components. For the 
two-particle final state one obtains 
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a!r\Born, loop) = ^ E 9^JHi^"'\Born, loop) (11) 



with a = U, L, I for i = 1 and a = F, A for i = 4 as above. 

Next we turn to the 0{as) tree graph contributions. The relevant three particle final 
state phase-space is given by 

2 2 

V q V q 

^^^3 = \2 dcos9dx g dydz —- dcosO dydz (12) 

where the transition to the last expression is again due to the azimuthal integration. We 
have introduced two phase-space variables y = 1 — 2pi ■ q/q^ and z = 1 — 2p2 ■ q/q^- The 
0{as) tree graph contributions to the various cross sections o"^™^ are written as 



111 4 

Tra f q 



3q'^ \ IGir^v ~[ 



J2g,, J Hi'^"^\y,z)dydz\ (13) 



with a = U, L, I for i = 1 and a = F,A for i = 4, as before. It is convenient to introduce 
the tree graph helicity structure functions Hl^"^\tree) by defining 

Hi^"^\tree) = -f^ f Hi^"^\y, z)dydz. (14) 

The Born term and the 0{as) corrections Hl^™'\Born) and Hi^'^\as) = Hl^"^\tree) + 
Hi^'^\loop) will be referred to as the unpolarized and polarized structure functions to 



^It is important to keep in mind that the transverse and normal spin components are defined w.r.t. the 
beam frame. When defined w.r.t. the event frame the transverse and normal spin components average to 
zero after azimuthal averaging. 
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leading (LO) and next-to-leading (NLO) order, respectively, while the sum of the LO and 
NLO contributions will be referred to as the 0{as) results. 

In summary, one has three unpolarized and seven polarized hadronic helicity structure 
functions. It is instructive to list them together including a specification of whether they 
are fed by the parity conserving {pc) or by the parity violating (pv) part of the product 
of hadronic currents and to which of the two classes of the so-called T-even and T-odd 



structure functions they belong to. One has 

unpolarized: Hu{pc), Hl{pc), Hpipv) T-even (15) 

longitudinally polarized: H[j{pv), H^{pv), Hp{pc) T-even (16) 

transversely polarized: H^{pc), Hj{pv) T-even (17) 

normal polarization: Hf{pc),H^{pv) T-odd (18) 



If one neglects contributions proportional to the imaginary part ImXz of the Breit-Wigner 
line shape of the Z-boson (see Appendix A) the T-odd helicity structure functions (pv) 
and Hf{pc) are contributed to by the imaginary parts of the one-loop amphtudes leading to 
nonvanishing triple product correlations of the type Sf{ly<pt), whereas the T-even structure 
functions obtain contributions from the Born term, the 0{as) tree graph contributions and 
the real part of the one-loop contributions. 

If one includes the contributions proportional to the imaginary part Imxz the structure 
functions Hp{pv), H^ijrv) and Hj{pv) are also contributed to by the imaginary parts of the 
one-loop contributions, and, vice versa, (pv) obtains also contributions from the Born 
term, the 0{as) tree graph contributions and the real part of the one- loop contributions. 
All the latter contributions originate from the {VA — AV) part of the product of hadron 
currents and thus belong to the class of helicity structure functions if^^™) according to 
the classification of Eq. ([2]). The latter contributions can only be probed through the 
imaginary part of the Breit-Wigner resonance shape which is strongly suppressed for {tt) 
production. In fact, the contributions coming from the imaginary part of the Breit-Wigner 
resonance shape are of order 0{lmxz{.<f) / ^^Xzi^f)) and can thus safely be neglected 
for top quark pair production. For example, in the threshold region of top quark pair 
production Imxz/ R^Xz is approximately 0.1% and decreases further with a 1/g^ power 
fall-off behaviour. We shall nevertheless include all contributions for completeness 

and for possible applications in (bb) production where the contributions cannot be 

neglected in the Z resonance region. 

3 Covariant expressions for the projectors 

The projectors Ua will be written in covariant form. We go to the rest frame of the gauge 
boson such that q = {\/q^; 0, 0, 0). The z axis is defined by the momentum direction of the 
top quark. For the top quark momentum one has 

Pi = ^ (i - 0, 0, V(i-y)'-e) (19) 
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{y = for two-body decays) with ^ = 1 — f ^ = Am? /q^. We construct a four-transverse 
quark momentum and a four-transverse metric tensor 

gfMu = gf,u - PiM = 9f^uPi = Pif, ^Qf, (20) 

and use q and pi to build up two elements of a coordinate basis, 

Q.''' / y(f ) ( = (1) 0) O5 0) the gauge boson rest system ) , (21) 



63 = \f>i/ \l {pi ■ qY/q"^ — ^ = (0; 0, 0, 1) in the gauge boson rest systemj . (22) 

In covariant form the longitudinal spin vector of the top quark reads (see e.g. [15]) 

s'^ = -{q" - ^p',)/^{piq)ym^ - q^. (23) 
m"^ V 

In the gauge boson rest system Eq. ( l23l) turns into 

/ = i= (v/(i-y)^-e; 0, 0, 1 - , (24) 

while in the top quark rest system one has = (0; 0,0,1). The longitudinal spin vector 
can be seen to be a linear combination of the two basis vectors cq and 63 and does not 
provide a new direction in our vierbein basis. The projectors that can be constructed with 
the help of Cq and 63 are limited to the three projectors 

np = i£ per (^3(^0 (25) 

where e^^p^j is the totally antisymmetric Levi-Civita tensor with £0123 = 1- They project out 
the three unpolarized and three longitudinally polarized helicity structure functions where, 
according to Eq. ([3]), the polarized structure functions H^^ (a = U, L, F) are obtained from 
iff = i7r {Hl,is')-Hl,i-s')). 

The transverse and normal polarization vectors of the top quark are defined in the 
beam frame. Viewed from the event frame they are given by 

ct = (O;cosx, -sinx, 0), cat = (0; sin x, cos x, 0). (26) 

These two vectors therefore allow one to span the beam plane and a plane perpendicular 
to the beam plane in event frame coordinates. With these new elements it is possible to 
construct the remaining additional projectors. They read (m = T, A^) 



nA{^m) — -^^-^^p-vpuG^oKni ^ {em) — '^^^'^m^S ~ ^S'^m)- (27) 



For example, according to Eq. ([3]), one obtains the structure function Hf^ by calculating 
Hf = n>^''{eT)(H%{s^) - Hl,{-s^)). Note that since i7,(e^) = iTKer), i7/(eT) = 
— i7|(e7v), nA^e^) = — i7^(er) and i7yi(er) = i7^(eAr) the primed projectors are redundant. 
This set of four (Eq. fl27|) ) and six (Eq. fl25|) ) covariant projectors allows one to calculate 
the complete set of ten helicity structure functions from the hadron tensor. 

In the following we list the Born term and loop contributions calculated already in 
previous papers [HI M, [IQ]. The nonvanishing unpolarized Born term contributions are 
given by 

Hlj{Born) = 2N^q^{l + v^), Hl{Born) = N^q^{l - v"^) = Hl{Born), 
H^{Born) = 2N^q^(l - t;^), H^{Born) = AN^^v. (28) 

The longitudinally polarized contributions read 

Hfj\Born) = AN,q%, H],\Born) = 2N^q^{l + t;^), 

Hl\Born) = 0, Hf{Born) = 2N,q\l - v^). (29) 

For the transverse and normal polarization components one has [9] 

Hf{Born) = N^q^^j^, H J {Born) = N^q^^j | = HY{Born), 




HY{Born) = N,q\^l^. (30) 

Note that one has Hi = Hi H}j = H}f, Hfj = Hf, Hj. = H^f, Hf = Hf, and 
Hf = Hf at the Born term level. We will return to these relations when we discuss the 
0[as) tree graph contributions. 

Note that the transverse and normal spin components T and N are proportional to 
f^ = 2m/f(p. The origin of this suppression factor is a helicity flip suppression factor at 
the j/Z — ti vertex. The same suppression factor also occurs in the 0{as) one- loop and tree 
graph radiative corrections to be treated later on. It is clear that this overall suppression 
factor is not important for (tt) production in the threshold region and not very significant 
in the range of beam energies considered in this paper. Altogether this means that the 
transverse and normal spin components of the top quark are non-negligible in the present 
application P, |TT] . 

Most of the nonvanishing one-loop contributions have already been given in [HI [9l [TO] 



Hlj[loop) = 4N,q^ ((1 + v^)ReA - 2v'^ReB) , 

Hliloop) = 4N,q^ ((1 - v^)ReA + 2v^ReB) , 

Hliloop) = 2N,q^ [{l-v^)ReA + v^ReB) = Hl{loop), 

Hl{loop) = -SNcq^vImB, 
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AN^q^ ((1 - v'^)ReA + 2v'^ReB) , 
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Ncq^vJ I {2ReA + (1 - 3^)ReB/^) 
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N,q'vJ^il+OlmB/C 



(31) 



where the real part of the form factor A and the real and imaginary parts of the form 
factor B read 
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(32) 



The imaginary contributions Hp{loop) and H^{loop) complete the list of one-loop contri- 
butions given in [HI [HI [ID] • We are now in full agreement with the one-loop contributions 
given in [11]. ImB contributes to the T-odd structure functions Hf^ and as men- 
tioned after Eq. f|T8|) . The infrared singularity has been regularized by the introduction 
of a gluon mass ma via = Ag^. The loop induced infrared singularities in the real 
part of the one-loop contributions can be seen to cancel against the corresponding infrared 
singularities in the tree graph contributions to be treated later on. 

In the next section we will present our results on the cut-off dependent helicity structure 
functions. They must coincide with the fully integrated results written down in Refs. [H 
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[9l [To] when the cut-off is taken to its maximal value £'G(niax) = (g^ — 4m^)/(2v^)- This 
will be verified in Sec. 5. 



4 Exact result up to a given gluon energy cut 



In this section we present the results 
of our calculations for the 0{as) cor- 
rections to the helicity structure func- 
tions with a given cut on the gluon en- 
ergy. We define a scaled gluon energy 
cut A = Ec/y/q^ and do the phase- 
space integration in the region < 
Eg < Aniax"\/?- The maximal value z 
that the cut parameter A can take is 
Amax = (l~0/2- In terms of our phase- 
space variables y and z the cut phase- 
space is defined by < y + z < 2A. 
In Fig. 2 we have drawn a {y,z) phase- 
space plot choosing a specific value for 
^ = 0.1 for illustrative purposes. The 
shaded area corresponds to the integra- 
tion region with the specific choice of 
cut value A = 0.3. The upper bound- 
ary of the integration region is given 
by the straight line z = —y + 2 A. 




0.3 0.4 



.8 0.9 1 



y 



Figure 2: Phase-space diagram with gluon cut 



The full phase-space is bounded from above and below by the two functions 2;+ and Z- 
where 



2y-2y^-iy±2yJ{l-yY-i 



The upper gluon cut given hy z = —y + 2A intersects the two boundary curves fl33|) at 



(33) 



Z/i = A 1 - 



'1-2A-^ 
1 - 2A 



2/2 = A 1 + 



'1-2A-^ 
1 - 2A 



(34) 



Since the phase-space is symmetric with respect to reflections along the diagonal, the 
corresponding values are zi = 7/2 and Z2 = yi- 

From a visual inspection of the phase-space plot Fig. 2 one can see that one has to 
discuss two cases when integrating the cut phase-space depending on whether (case A) 
A < Atrans OT (case B) A > Atrans- The transition value Atrans = (1 ~ v^)/(2 — v^) is defined 
by the A value at which the straight boundary line of the cut intersects the phase-space 
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boundary at the point (?/2, -22) = (1 — V^, a/^ (1 — V^)/(2 — V^))- At this point the tangent 
of the full phase-space boundary is vertical. From an inspection of the phase-space plot 
Fig. 2 one concludes that in case A the integration region is divided into two parts, whereas 
one has to consider three integration regions in case B. The specific example shown in Fig. 2 
corresponds to case A. 

Let us denote the general y- and z-dependent tree graph integrand in case A by I{y, z). 
One has to do the two integrations 

yi rz+ ry2 r'2X-y 

/ Iiy,z)dydz+ / / I{y,z)dydz, (35) 

while in case B one has an additional integration, viz. 

ryi ry2 r2X-y f^^V^ 

/ / I{y,z)dydz+ / / I{y,z)dydz+ / / I{y,z)dydz. (36) 

Jo Jz— Jyi Jz- Jy2 J z— 

It is clear that one should recover the fully integrated results listed in Sec. 5 when setting 
A to its maximal value Amax = (1 0/2- When comparing to the fully integrated result 
one has to discuss case B with A = A^ax = l/i = I/2 = (1 — 0/2- In this case the second 
integral in Eq. (|36i) vanishes and the remaining two integrals can be merged to give 

I{y, z)dydz (37) 



which corresponds to the fully integrated tree graph contribution entering the full NLO 
result given in Sec. 5. 

Let us return to case A involving the two integrations in Eq. fl35l) . For most practical 
applications case A will be the relevant case since the ratio 

'^trans 2 (QQ\ 

remains close to 1 over most of the range of ^ values. The integration over z is straight- 
forward. The second integration over y is done by using the Euler substitution 

y=l-i(\±4, (39) 

Eq. fl5I?]) is easily inverted. The y-integration limits y = 0, yi, y2, 2A, and 1 — a/^ translate 
into w = Wo, wi, W2, w\, and 0, where 



Wo 



1/1,2 - 



Wx 



The value wx corresponds to the intersection of the upper gluon cut boundary with any 
of the two axes. In addition to the velocity parameter v = y/1 — C, we introduce modified 
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velocity parameters Vi — ^(1 — yj)^ — ^ and vx — ^J{1 — 2A)^ — ^. We shall also use 
the abbreviations a — 2 + y/^ and b — 2 — y/^. Our results for case A read {N — 
asN,Cpqy{4nv)) 



Hi, = iv|2(2 - 0'\to- - to+) - (8 - lOe - a(^i- - ^1+) + 

+ ^(1 - V^)(2 + - - (6 - llO^' + 2(2 - 04+ - + 

-^^o , lot t2^/) , 0\2 , t2 



- -(8 + 12^ - e')4+ + I 16A - 8A^ -2^ + e- ) ^7+ + 



^16X - 6A2 - 2(1 + 4A + A2)e + - ^) 4+ + 2yi(4 - |/i)4 - 2^2(4 - ^2)4 + 
2(1 - 2A - (1 - X)V^) 



(1 - v^)v^ 



(l-2A)(2-e) + 



+ 2A^ - (4 + 3A)e^ + 3e - rr^^) ^3 + 



4^24 + 5CK + ^^^-l(48 + 50.. + 6.? + 

^ ^Ar|2(2 - 0(^0- - ^0+) - (4 - Oih- - h+) + ^fa^ - + 2^ + 2^4+ + 



- + (8A - 0(4+ - 4+) + -^4+ + 4z/i£i - 4^24 + 

- 2 ^~^^^J^j^^^ (l - 2A - e + A^)4 - 5^;i - 5y, - 5v, + 5y,^ 

< = 7v|4(2 - ^)v{to-. + to+) - (8 + 2e + )(ti- + ti+) - ^(1 + ^)(2 - + 
+ 4v£4- - 3(2 - 50^;4- + ^(28 - 1704- + I (6 - 16A - 130^^ + — J 4- + 
-Vx f6-4A-13e-2Ae + ^(7A + 2A2 + + ^) 4- + 



(42) 
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_(l_6A + 8A^-9e + 3AOv^)(^-^) + 



+ \; -2V^^''^ { ^ - «A + 8A- - 7^ + 5A^ + 2A-^ - 3^' + 



+ (1 - 6A + 8A' - 9^ + 3A0v^|k - «;2) + 



4^ ' ^ y^J 2{b-awi) 2{b + aw2) 

+ ^(24 - SSe - 8v,)y, - 2yf + ^(24 - 33^ + Sv2)y2 - 2yi^ (43) 

= Ar|e(2 - 0(to- - to+) - 2^(2 + Oih- - - v^(l - ^)(2 + 4^^ - + 

+ 1(16 - 54^ + Se)v + Ch+ + 16^4+ + 4(8 + 8^ - 3^')4+ + 
4 Id 

- -^(8A - 28A2 + 16A3 + 16A^ + ^ - 12AC - 8X^^ - C^)^^^ + 

- ^(8A + 4A2 - e - 8Ae + - e (21/1 + ^yi) 4 + i (22/2 + ^2/2) 4 + 

+ 2A - ;eKr-'!l^ ^^ - 8A + 8A- - e + 2Ae - 4A-e + 7,^ - + 
- (2 - 6A + AX' + Se - 3Ae - 2A2e + 3^)7^)4 + 

- 1(32 - 72^ + 5^^ - SMv. + - 1(32 - 72^ + 5^^),, + 

= eA^|(2 - 0(^0- - io+) - 2^^(ii- - - ^(1 - ^fOtw + ^(22 - 30^ + 



(44) 



1 1 2X'^f 

1+ + -(8 - 8e + 3e%+ + ::{SX + 4A^ - 0(4+ - 4+) - ^4+ + 
ib 2 V 

22/1 + ^2/0 4 - (22/2 + ^2/2^) 4 + 2 ^~,^,^_~^;^r^^^ (l - 2A - C + AV^)4 + 



(i-ve)ve 

-1(72-5^ + 8.0.1-^^/^-1(72-502/1 
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e(10 + 30(^1- + - |(24 - 704- - 13<(4- - 4-) - ^4- + 

^1 2(l-2A+v^) 2(1-2A-Ve) J' 

I 3(l + v^)^(2 + v/C)v^ 3(l-v^)^(2-v/C)v^ V, , 
"'H^^ 2(1-,, + 2(l-,,-Ve) )''^ 

+ r+ 2(l-,. + V!) 2(l-,.-V^) J'^ + 

+ 1 2 - 8A + 8A2 - 7e + 5Ae + 2A2e - 3^' + 



(1-2A-V?)V? 



(l-6A + 8A2-9e + 3AOye| 



Wi W2J 



_ l^^±iL^l^/ 2 - 8A + 8A^ - 7e + 5Ae + 2A^^ - 3^2 + 
(l-2A + Ve)Ve I ^ 

+ (1 - 6A + 8A' - 9^ + 3Ae)v^|(«;i - W2) + 
+ (2 + 70^1 + (2 + 70yi - (2 + 70^;2 + (2 + 70^2) (46) 



Hi' = iV|2(2 - iY{to- - to+) - (8 + 2e + - h+) + 

- 2(6 + 0^ + 2(2 - 04+ + 8^4+ - ^(4 + 6e - 3e')4+ + 
+ (l6A(l - A^) - 2(1 + 6A)e + e- 4+ + 

- (^2A(8 - 3A) - 2(1 + 4A + A^)^ + - ^ j 4+ + 

+ 4(1 - 2A)(3 - 2A)4+ + 2yi(4 - 3^ - yi)4 - 2^2(4 - 3^ - ^2)4 + 

- (15 - 22A + 8A' + e + A0\/e)4 + 
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_ 3(4 - 7£ - 4,.).. - + i(16 - 21?),. - d] (47) 

Hj.' = ejv|2(2 - 0(i„_ - i„+) - (4 - 0(f,_ - + ^{1 - sfQt„ + 2v + 2U^ + 

- + (8A - 0(^7+ - 4+) + ^4+ + - 4^/24 + 

- 2i:i^^ii^-^(l - 2A - e + Aye)4 - 5^;i - 5|/i - ^v<, + 5|/2| (48) 

= 7V|4^;(2 - 0(to- + ^0+) - 2(4 - 50 (ti- + - ^(1 + Ve)(2 - ^C)' + 
+ 4^;£4_ - 6v4_ - 8C4_ - — (3 - 14A + 20A2 - SA^ - 2^ - + 

+ ^(3 - 8A + 4A2 - 2^ + SAC - + 4(1 - 2A)(3 - 2A)£9- + 

- 2(3 - + 2(3 - ^2)^^24 + 

. 1(12. 16A.,-4.0..^-^. 1(24.0.-..^. 

- 1(12 + 16A + e - 4.2).. + + ^ + + Oy. - y.^} (49) 

- 1(1 + Ve)(2 - Ve)' + 2^;£4- + 1(72 - 30e - 3^)4- + 



+ ^;,(|-4(l-2A-0(l-4A-0 



3(l + v/e)^(2 + Ve) , 3(l-v/e)^(2-Ve) \. , 
2(1-2A + Ve) 2(1-2A-V^) J'-^ 

- 1(12 + 70^(4- - 4-) - 1(8A - 4A^ - e - SAC - A^^ + f )4- + 2(1 - 2A)4- + 

(i , 3(1 + v/!)^(2 + 3(l-v/e)^(2-vT) \ 
+ I2+ 2(l-,, + Ve) + 2(l-,,-Ve) Y'^ 

(i^ 3(l+v^)^(2+v^) 3(1- v^)^(2-v/C) \ , ^ 
^^2+ 2(l-y2 + Vf) 2(l-y2-Ve) T 
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+ ^ J^f [ 2 - 8A + 8A^ - 7e + 5Ae + 2A^e - 3^^ + 



_ (1 _ 6A + 8A^ - 9e + 3X0 (—-—) + 

^ J \Wi W2J 



+ ^ (^^^ J^^^f [ 2 - 8A + 8A^ - 7e + 5Ae + 2A^e - 3f + 

+ (1 - 6A + 8A' - 9^ + 3A0v^|k - «;2) + 

- -(28 + 50.1 - ^ + - g (28 + 50yi + 

+ 1(28 + 50.. + ^ + ^^^^ - ^(28 + 50^4 (50) 
8 y2 4(0 + aw2j 8 J 

= ^ Vf ^{2(2 - 0(^0- - W) - ^(16 - 30(ti- - + 

- ^(1 -fm- im + ii)t^ - ^(16 - 30. + 

+ 2^4+ + 2(7 - 04+ - ^(8 - + 3e')4+ + 2(1 - 2A)4+ + |yi4 - \y2^2 + 

- e (1 - A - 4+ - (|4A(2 - A) - e + ^) 4+ + 

Hf = i/|A'{2(2 - Q(i„^ - W) - i(8 - 30(«i^ - *,+) - 5(1 - sfOU, + 
+ 5(4 - + 24+ + 2(1 + + i(2 - 0(4 - 35)4+ + (8A - 5 - A5)4+ + 

- [i\(2 + A) - 5 + 4+ - 2(1 - 2A)4+ + 

+ i(8 - QsiCi - i(8 - 0»24 + ^ " ^V- ^ ^''^ '^ - 2A + Ay/i - 04 + 

-5(^^-««'"-4(F=^-5(«2-««^' + 
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+ 1(2 - V^)2(l + v/e) + 2t;£4- - ^(8 + 0^4- + ^(8 - 30e + 3e')4- + 
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- 2(1 - 2A)£9- - ^(8 + 0^1^ + ^(8 + 0^2^ + 

+ 1(52 + 50.. - ^ - + 1(52 + 50.. + 

8 Hi 4(o — awi) 8 

- 1(52 + 50.. + ^ - j^^^^ + 1(52 + 50.4 ^^3) 
8 .2 4(6 + atus) 8 J 

The logarithmic rate terms ii and the double and dilogarithmic rate terms to-? ^o+) ti-y 

and are listed in Appendix C. Note the exact 0{as) tree graph relation Hp^ = 
which was also noticed in [H]. We have not been able to derive this relation from general 
principles. 

We shall not dwell on the technical details of how the finite integrals have been cal- 
culated but rather concentrate on the class of IR-divergent integrals. For instance, the 
integral 

/.(-I, -1) = r r ^ = r in f ^ (54) 

^ ' ' Jo Jz^ yz Jo \z-{y)l y ^ ' 

is IR-divergent and will be regularized by a gluon mass mc = y/Aq'^. The introduction of 
a gluon mass changes the lower y limit from to = A + \^7i^, and the z limits to 

1 



= - - + + 2^ ± - A)' - (1 -y?-^]- (55) 

Therefore, the integration over z gives rise to 

m , (2y-2y^-iy + 2Ay + 2A + 2J {y - Af - A^J {I - yf - i\ dy 
-1 = / In ' ' 



'y- \2y- 2y2 -^y + 2Ay + 2A- 2^{y - A)2 - Aev'(l - y^ -^J V 

(56) 

This integral is not analytically calculable for general values of A. However, we can divide 
the integral into a divergent and a convergent part which are separately calculable as 
long as A is a small parameter. The residue of the divergent part should coincide with the 
residue of the original integrand at the IR-singular pole at . = 0. A simplified IR-divergent 
part can be constructed from the full integrand by neglecting higher powers in y whenever 
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possible. Before this approximation we shift the integration by —A in order to facihtate 
the expansion around the lower boundary. We obtain 



/?(-!. -1) = r f i!+^"^+f-%I5 1 (57) 



/Xi \ (1 + v^)?/ — 2f a/?/2 — A^y y 
This integral can be calculated analytically and one obtains 



In the case A ^ we have the limiting value (we write e = v^A^) 



/f(-l,-l)^21n(i±ii)li,„r* (59) 



1 — vj e^OJe y 

which is an ill-defined quantity for e = 0. However, we can subtract the singular piece 
from the original integral also taken in the limit A ^ 0. As a result the divergences cancel 
and one obtains 



yi 



l,-l) = lim<!/ In I \ = I — - 21n 



2-2y-^ + 2J{l-yy dy ^, fl + v\fyidy 



s^o]^Je y2-2y-i-2^{l~yy ~ij V \l-vJJe y 

(60) 

or symbolically 

(-1, -1) = Jim {/r (-1) - 2 In ([±^) /^(-l)} (61) 

where the primes indicates that the lower limit is e. With the Euler substitution Eq. ( l39i) . 
and after partial fractioning according to 

dy dw dw dw dw , , 

+ + 1 q (62) 



y Wo — w Wo + w 1 — w 1 + w 



one obtains 



rba/ 1 



C'(-l) = In 



2-2|/-e + 2J(l-|/)2-A dy 



2-2y-i-2^{\-yf-ij V 

= lo-iO - itM - loliO + lolM + 

-Itiw'o) + Itiwi) + I'MO - (63) 
^^(-1) = [ — = lo-iw'o) - Io~iwi) - lo+iw'o) + lo+iwi) + 



y 



-h-{w'o) + h-{wi) + h+{w'o) - h+{wi) (64) 
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where 



^ uJl-i4 + -..V (65) 



The variable w'q tends to wq for e ^ 0. It is instructive to note that the divergences now 
reside in the terms I^tiw^) and Iq^{w'q) which contain the integrand factor {wq — w)~^. 
We obtain 

= t'^{w) - 2 In ([^) Inlu'o - w), h.{w) = - \n{w^ - w) (66) 

where tj, is a decay rate term which vanishes in the hmit w ^ Wq. For this reason the two 
expressions in Eq. (|66|1 do not contribute to the convergent part at alL Using Eqs. fl63|) 
and flMl) we can calculate the convergent part in Eq. fl6T|) and add the divergent part in 
Eq. ([58]) to obtain 



/,(-!, -1) = tp-ln(^l±^)lnA-|/o^^l(«;i)-21n([^)/o-(t^i)} + ... = 

= C-(^o) - to-K) - C+K) + + • • • - In (^) In A. (67) 

The dots indicate further contributions according to Eqs. (163!) and fl64l) where we can 
replace Wq by Wq. The decay rate terms t^Q_{w), t^^j^iw), . . . are listed in Appendix C. It 
is obvious that to-('^o) in Eq. (1^ has to be replaced by the special value tp defined in 
Eq. dSH]). 

We now turn to case B when A > Atrans- As discussed in the beginning of this section 
this entails the calculation of the second integral in Eq. fl36l) which has to be added to the 
first and the third integral in Eq. fl3^ . The latter two integrals are already known from 
case A. Using some additional decay rate terms listed in Appendix C the results for this 
additional phase-space portion are given by 

Hh = Ni^2{2-0\t^o--to+)-i8-m-eKtL-tl^) + 

+ \/e(l - v/?)(2 + - 304 - 8(2 - 0<+ - IQ^^U + 

+ (1(8 + 12e - e) + 82/2 - 2ytl + ^(24 + 50^2 - f^'f^ }' ^^^^ 



Hi, = eAr|2(2-0(to--iS+)-(4-0(tt-t^+) + 

+ V^(l - V^)4 - 8<+ + + %2) ^2 + WV2^, 
Hfj' = Ar|4(2 - Ovi^o^ + - (8 + 2^ + 3e')(t^_ + t^+) - 8(1 - 0^- + 

+ 16(1 - 0^1- - ^(28 - I70mi- + ^eJ - 16(1 - O^r^ + 



(69) 
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+ ^(64 + 64v^ - 352^ + 232^^^ + e) - ^(8 - IIO2/2 + ^vl - |^^}' (^0) 

Hi = 7v|e(2 - 0(iS- - - 2^(2 + 0(iL - ^^+) - ^(2 + 2^ - 7^ + 3^^)^ + 

- 4e<+ + 16e£^+ - ^ (^1(8 + 8^ - 3^') + 2^2 + \yt) ^2 + 

+ i(32 - 72e + 5e^-8e.>.-^^), (71) 
Hi = eiv|(2 - 0(^0- - ^0+) - 2(1 - CM- - - V^(l - V^)C - 4<+ + 

- (^(8 - Se + 3f ) - 2,. - Ivl) + i(72 - 5e + 8,.).. + |^), (72) 
Hf = 7v|e(10 + 3O(i5_ + t?+)+e(24-7O^L + ^e(24-7O^^_ + 

+ (^(2 -im- iif - ^(2 + im + Ve)^ - 4^.2) + 

+ 2(1 - ^)(2 - + 130 - 2(2 + 7e)y2|, (73) 

= iv|2(2-0'(t^_-t^+)-(8 + 2^ + a(^?--^5+) + 
+ ^(8 - 10^ + e + - 8(2 - 0<+ + 8^^^+ + 

+ (^(4 + - ?>e) + 2(4 - 3C)|/2 - 2yi) l\ + ^(4 - 7^ - ^y<^)v^ - |^^^}' (74) 

E^^ = eiv|2(2 - 0(t^_ - tS+) - (4 - 0(tL - + 

+ \/e(l - iiYl - 8^^4+ + + 41/2) + 10t;2|, (75) 

E% = Ar|4(2 - 0^(iS_ + 2(4 -50(i5_ + t5+)- 8(1-0^^- + 

+ m%_ + 8e^^_ - 16(1 - - 2(3 - 1/2)^2^2 + 

+ ^(80 - 64^ - 8e + e) - ^(24 + 0^2 + 2^2^ - (76) 
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- i(4 - 0(10 + mi- - ^(72 - 30^ - + 

- 8(1 - m_ + + ^(2 - v/^)(l - fif + ^(2 + V^)(l + yj)^) + 

1 r r 1 ff4-f)^ 1 

- -(304 - 496 + 208^ - 2A^i^ - f) + ^(28 + 50z/2 - (^7) 



HY = ^y|iv|2(2-0(tS--io+)-^(16-30(t?--t?+) + 

- 1(8 - loy^ + e + - 8<+ + 2(7 - + 

+ Q(8 - 6e + 3e^) + le..) + ^(4 + 50.. - l^j, (7^ 
hT = ^yf ^{2(2 - 0(^0- - to+) - ^(8 - SOitl_ - tl^) - - V^)4 - 8<+ + 
+ 2(1 + 0^5+ - f ^(2 - 0(4 - SO - 1(8 - O1/2) + 7(52 - 50t;2 + " ^^""^ ^ 



^v"^ "s;^. 4(42/2 + 01' 



(79) 



HT = ^^iv|2(2-Ot;(tS_+tS+)-^(8-130(t5_+t^+)-4(l-O^L + 
+ ^(24 - 2^ - 3e)it - 1(8 - 30^ + 3en- - 8(1 - 0^7- - ^(8 + 0^2^2 + 
+ 1(208 - 208V^ + lee - 24^ ' " l^^^ + ^Ol/. + ^^}- 



^0) 



Note again that one has the remarkable relation H]f = H^. We mention that, differing 
from Ref. [12], we have been able to obtain a closed form result for the cut-dependent 
structure function Hp (see Eqs. (H9|) and (!76|) ). 

Numerically, the contribution of the second integral in Eq. (!36|) calculated above is quite 
small. This is because the relevant integration region is far away from the IR region where 
the rate is largest. Nevertheless, this contribution is needed if one wants to check on the 
consistency of our case B result with the fully integrated results in Refs. [3 El El HOI [11] . 
In fact, we have performed an explicit check that for each of the unpolarized and polarized 
rate functions the sum of the three integrals in Eq. (I36l) reproduce the full phase-space 
result calculated previously in Refs. [3 El O jTUl IH] when the gluon energy cut is set to its 
maximal value Amax = (1 ~ 0/2 (which corresponds to setting yi and 2/2 to (1 — 0/2 in 
Eq. We have also checked that our exact result converges to the soft-gluon expression 

to be derived in Sec. 6 when A 0. 
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5 Fully integrated 0{as) results 

The cut-off dependent helicity structure functions calculated in the previous section must 
coincide with the fully integrated results written down in Refs. [HI El ID] when the cut-off is 
taken to its maximal value. For the convenience of the reader we collect the fully integrated 
results of [U [9l [10] and list them in terms of the sum 
As before we define (A^ = asNcCpq'^ /{4:Ttv)). One has 

Hhias) = Ar|(2 + 70t; + ^(48 -48e + 7at3 + V^(2-70t4 + 

+ e(2 + 30(t4 - h) - 2(2 - ((2 - 0(^8 - h) + 2v{ho + 2ti2)) }, 

- 2 ((2 - 0(^8 - tg) + 2i;(tio + 2t,2)) }, 

Hlias) = iv|^(16 - 46^ + 3^')^; + ^(88 - 32^ + 3e')t3 - ^^(2 - 70U + 

- e(2 + 30(t4 - ts) - e ((2 - 0(^8 - tg) + 2^7(tio + 2^12)) }, 
Hlias) = eiv|^(10-Ot^ + ^(24-16e-3ai3-v/et4 + 

+ ^{U - h) - ((2 - 0(^8 - tg) + 2v{tio + 2tu)) }, 

HUc^s) = iv|-8v^(l-V^)-8(ti-t2)+ 4(2-3^^3 + 

- 2(4 - 50^6 - 4^; ((2 - 0(^8 - h) + 2v{t,o + tn)) }, (81) 

Hlfias) = - (2 + 350 + v/e(8 + 290 -^(32 - 60^ + 17^(^1-^2) + 

+ 2(4 + 90t;t3 - (8 + 2e + 3e)t6 - 4:v ((2 - 0(^8 - ^7) + ^^(^10 + ^ii)) }, 
Hl\as) = 0, 

Hfias) = Ar|2(2 + 190 -27^(8 + 130 + 

- ^024 - 70(ti - h) - 2Q^vt, + + 30*6}, 
H'Aa,) = Ar|-2(2 + 30t^ + (24-12e + a*3+y^(8 + 0*4 + 

- - 0(^4 - ^5) - 2(2 - ((2 - 0(^8 - h) + 2v{t,o + 2tu)) }, 
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- 2 ((2 - 0(^8 - ^9) + Mho + 2^12)) |, (82) 



HfM = -^^^iv{48 + 17^-^(62 + 30 -^(4 -0(10 + 30(^1 -t2) + 

- 2(21 + 20^;t3 + (16 + 70^6 + ((2 - 0(^8 - + 2v{ho + ^ii)) }, 

HY{as) = -^yfiv{(8-30t'-^(72- 38{ + 3at3- V^(10-O^4 + 
+ (8 + 0(^4 - is) + 4 ((2 - 0(^8 - h) + + 2ti2)) |, 

HIM = -^yfA^{-(20 - 30^-^(32-l4e-3a^3-e^/e^4 + 

+ i{U - ts) + 4 ((2 - 0(^8 - h) + 2^;(tio + 2^12)) }, (83) 



2 V 2 



= ly|Ar|20 + 9e - ye(26 + 30 - J(24 - 2e - 3e^)(ti - t2) + 

+ 2(1 - 60^;i3 - (8 - 130*6 - ((2 - 0(^8 - ^7) + 2v{tio + hi)) }, 
HTM = l\llNil + Ovn. (84) 

The fully integrated 0(0;^) results are given in terms of the rate functions ti to ti2 which 
are listed in Appendix B. It is clear that one again has the relation H'p[as) = Hij[as) 
because both loop and tree contributions satisfy this identity. 



6 The soft-gluon approximation 

The basic ingredient of the soft-gluon approximation (SGA) for the tree graph matrix 
elements is the eikonal approximation where the gluon momentum is neglected in the 
numerators of Feynman diagram contributions. In the eikonal approximation the hadron 
tensor is proportional to the Born term. In the present case one has 

HUsoft) - i.e. - . j£^) H^Born, (85) 
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where W^^{Born) refers to the Born term tensor in the two-body case where q = pi + P2- 
On the other hand, the eikonal factor multiplying W^^{Born) refers to the three-body 
case where q = pi + P2 + P'i and depends on the dimensionless three-body phase-space 
variables x = Ec/Vq^ = Psq/q"^ and u = {pi — P2)q/q^- When integrating HJ^^^soft) over 
the three-body phase-space the Born term contribution W^^{Born) can be taken outside of 
the integral. In this sense the integration on the soft-gluon factor in Eq. fl85|) is universal 
in the sense that it is process and polarization independent. 

When projecting the eikonal contribution in Eq. fl85|l onto the various helicity structure 
functions one recovers the various Born term contributions Hl{Born) listed in Sec. 3. 
Referring to the integration measure in Eq. ( |T^ and using dy dz = 2dx du one obtains 

where 

^ (l-2x + A)(.^-(.-Ay^)+^(x-A)^ 
The limits of the u integration are given by ±u+ where 



1/2 

I — —I— /\ — r \ 

u+{x) = (x - A)- 



2 ^ , 1 - 2s + A - ^ 



2a; + A 



After integration over u one obtains 



Further integrating over the scaled gluon energy x from A to A one finally has 



+2i; (ln(— ) +21n( ^ ) ) -lnzo + 4AlnzA + 



ZqJ \ Z\Zq — 1 ^ 

1,__2/^aA . ^ (zxZq-1\ ,1, 2 



+Li2 + Li2 Li2(l - zxzo) - Li2(l - zl)\ I (90) 



where 



l+vj V Zq 



l + v yr^2A + v/l-2A-e 

^0 = ^ , zx = . 91 

l-v vl - 2A - vl - 2A - ,^ 



The function /ipii, will be referred to as the eikonal form of the SGA factor. 
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For A one obtains 



Following the literature p3l [T6j we shall refer to the SGA factor (1921) as the soft-gluon 
approximation of Eq. (l90ll . 

In addition to the check on our case A results discussed in Sec. 4 we have performed a 
second and independent check by taking the A — > limit in the relevant exact expressions 
in Sec. 4. In this limit the exact result can be seen to factor into a Born term contribution 
times the soft-gluon factor given in Eq. fl^^ . This proves that the exact results given in 
Sec. 4 have the correct soft-gluon limiting behaviour. 

In order to be able to compare the eikonal SGA factor Eq. (1901) and its approximate 
version Eq. (p2|) we (minimally) subtract the IR-divergent piece hj^ from both expressions 
where 

.„^-^{(2..(2-0.(i^))>„(-L)j. ,03, 

The remaining IR- finite pieces are then /igj^ = /leik — hm and /^sga = ^sga — hiR. In 
Fig. [3] we show a plot of the relative fraction {h'^^]^ — h'^Qx) I ^sga a function of the cut-off 
parameter A/Amax- Fig. [3] shows that \h!^^^^ > {h'^Q^l since both functions /i^ji^^ and /^sga 
are negative over the whole range of A. The SGA Eq. (p2|) is a poor approximation to the 
eikonal approximation Eq. (!90|) except for the region very close to the soft-gluon point. 
For ^/s = 1000 GeV the fractional deviation can become as large as 100% at the maximal 
cut value. 

As it turns out the eikonal approximation with the eikonal factor (l90l) approximates 
the exact result rather well numerically even up to the hard end of the gluon spectrum. In 
Fig. m we show a plot of the total rate a {=au+L) as a function of the cut-off parameter 
A/Ajnax for the three center-of-mass energies ^/s = 400, 500 and 1000 GeV where we take 
rrit = 175 GeV and as = 0.0964, 0.0941 and 0.0875, respectively, for the above three 
energies. The rates rise very quickly from the soft region to values close to the total rates 
showing that the contributions from the soft region dominate the total rates. The quality of 
the eikonal approximation becomes marginally weaker when the hard gluon region becomes 
larger with the increase of the center-of-mass energy. The exact result is hardly discernible 
from the eikonal result at the scale of the figure even for the highest cm. energy. The SGA 
approximation can be seen to be quite poor. Also shown are the respective LO Born term 
contributions which appear as dotted horizontal lines in Fig. |H The radiative corrections 
can be seen to be quite large. At the point where the 0(0;^) rate intersects the LO Born 
term rate the corrections go to zero. This can be seen to happen at A/Amax = 2 x 10~^, 
0.014, and 0.200 for the above three cm. energies. At even smaller cut values the total 
0{as) rate goes to zero altogether. This happens at A/Amax = 10~^^, 10~^, and 5 x 10"'' 
for the same three above cm. energies. It is clear that perturbation theory should not 
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Figure 3: Dependence of the relative fraction (/i^j^ — /isGA)/^SGA scaled gluon 

energy cut-off parameter A/Amax where Amax = (1 ~ 0/2- Curves are shown for the three 
center-of-mass energies ^/s = 400 (dotted), 500 (dashed), and 1000 GeV (full line). 

be used for such small values of A. This holds, in particular, for the polarization-type 
observables to be discussed later on since they are normalized to the total rate and are 
thus very sensitive to the vanishing of the total rate. It is important to keep in mind that 
the NLO rate goes to — cxd when A — > even if this is not apparent in Fig. |H 

In order to show the quality of the eikonal approximation in Fig. O we show a plot of 
the cut-off dependence of the relative difference of the exact cross section and the eikonal 
approximation (a — (Tcik)/cr for the same three center-of-mass energies. For ^/s = 400 GeV 
the relative difference is very small and remains below 0.1% over the whole gluon energy 
spectrum. For the largest energy shown {^/s = 1000 GeV), where the hard gluon region is 
the largest, the relative difference rises from zero at the soft end of the spectrum to about 
2% at the hard end of the spectrum. 

7 Numerical results 

Let us begin the numerical section by the statement that we shall, as in the previous 
section, always use a top quark mass of 175 GeV in our numerical results. Since all our 
results are given in analytical form the corresponding results for other values of the top 
quark mass can be readily calculated. For the strong coupling constant we take the same 
values as described at the end of the previous section. 

We shall divide our numerical results into two subsections according to whether the 
observables or structure functions have a nonvanishing or vanishing Born term contribution. 
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Figure 4: Dependence of the total rate (solid line: exact NLO; dash-dotted hne: eikonal, 
dashed line: SGA) on the scaled gluon energy cut-off parameter A/Amax where Amax = 
(f — 0/2. Also shown are the respective cut-off independent LO Born term contributions 
(horizontal dotted lines). Curves are shown for the three center-of-mass energies ^/s = 400, 
500, and 1000 GeV. 




Figure 5: Dependence of the relative difference of the exact cross section and the eikonal 
approximation on A/Amax where Amax = (l~0/2for center-of-mass energies y/s — 400 
(dotted line), 500 (dashed line), and 1000 GeV (full line). 
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Figure 6: Dependence of the ratio a {hard) /a (full) on A/Amax in the hard region where 
A denotes a lower cut-off. Curves are shown for the three center-of-mass energies y/s = 
400 GeV (dotted), 500 GeV (dashed), and 1000 GeV (full line). 

7.1 NLO corrections to nonvanishing LO observables 

We shall use a terminology where the NLO results are partitioned into a soft and a hard 
region by a cut-off value for the gluon energy Ec. The soft and hard regions are defined 
by their respective integration regions. In the soft region one integrates from zero gluon 
energy up to the gluon energy cut Ec including, of course, the one-loop results. In the hard 
region, one integrates from the (lower) gluon energy cut Ec to the maximal gluon energy 
E = {1 — ^)A/g^/2. We use this terminology to differentiate between choosing an upper 
cut-off (soft region) and a lower cut-off (hard region) even if the respective integrations 
extend into regions with maximal and minimal gluon energy. The hard gluon contribution 
can be obtained by subtraction. Thus, for example, a{hard) = a — a{soft). The definition 
of the two regions holds irrespective of the actual value of the cut-off energy. 

In Fig. [6] we show a plot of the ratio a (hard) /a (full) {a{full) = ex) as a function of the 
cut-off parameter A/Amax for the three cm. energies y/s = 400, 500 and 1000 GeV. Note 
that the hard gluon fraction is proportional to The hard gluon fraction is generally 
quite small. As the lower cut-off tends to zero a{hard) and thereby a {hard) / a {full) tends 
to -|-oo (due to the positive — log A singularity). Away from A = the hard gluon fraction 
then drops very quickly as the lower cut-off is raised and reaches zero at A/Amax = 1 
where there is no phase-space left. The hard gluon fraction becomes larger as the energy 
increases. For example, at A/Amax = 0.2 the hard gluon fraction is 1.5, 4.4, and 13.6% for 
\/s = 400, 500 and 1000 GeV, respectively. The corresponding soft-gluon fractions can be 
obtained by subtraction as mentioned above. 

We do not show corresponding plots for the other partial unpolarized and polarized 
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Figure 7: Dependence of the differential rate da /d cos 6 on cos^ in the soft region. Curves 
are shown for the three center-of-mass energies y/s = 400 (dotted hne), 500 (dashed hne), 
and 1000 GeV (full line) and three upper cut-off values A/A^ax = 0.2, 0.4 and 0.8 (from 
bottom to top). 

rates af™''' because they do not differ much from those shown in Fig. [61 This can be under- 
stood from the discussion in Sec. 6 where we demonstrated that the real gluon emission 
contributions are very well approximated by the eikonal approximation which in turn is pro- 
portional to the Born term contribution. This implies that all ratios af^\hard) / (j^^\full) 
are approximately equal to one another as well as approximately equal to a (hard) /a {full). 
An exception is af^ where the Born term contribution is zero. This case will be discussed 
in more detail later on. 

In Fig. [7] we show a plot of da /d cos 6 as a function of cos^ for the three cm. energies 
^/s = 400, 500 and 1000 GeV and for three respective cut-off parameter values of A/Amax = 
0.2, 0.4 and 0.8. The cos 6' dependence is marked and strongest for ^/s = 500 GeV showing 
that the forward-backward contribution ap is non-neglible. The radiative corrections are 
large for y/s = 400 GeV and y/s = 500 GeV similar to the total rate plotted in Fig. |H The 
cut-off dependence is generally quite weak showing that the bulk of the different partial 
rates comes from the region close to the soft-gluon point A = 0. 

In Fig. [8] we show a plot of ApB as a function of the upper cut-off A/Amax again for 
the three cm. energies ^/s = 400, 500 and 1000 GeV where we have defined the forward- 
backward asymmetry by 

^ a{forward) — a{backward) 
a (forward) + a (backward) 

Note that one has to separately integrate the numerator and denominator of Eq. ( IMl) 
over the gluon energy when calculating ApB- The radiative corrections are generally small 



30 



and the dependence on the cut-off A is quite wealc. ApB is largest for y/s = 1000 GeV 
as can also be appreciated by looking at Fig. [71 The radiative corrections are largest for 
y/s = 400 GeV. For example, for an upper cut-off of X/Xmax = 0.2 they amount to 2.7%. 

The radiative corrections to polarization-type observables p/™'' are in general quite 
small even if the radiative corrections to the polarized rates themselves are large. The 
reason is that polarization-type observables correspond to normalized density matrix ele- 
ments defined by the ratio of a polarized rate and the total rate. The radiative corrections 
to the numerator and the denominator tend to go in the same direction and thus tend to 
cancel out in the ratio. Take, for example, a generic polarization observable Pj;"^'' which, 
at 0{as), is defined bjifl 



P/")(0(a.);A) 



at\Born) + at\c^s;X) 
a{Born) + cr(as; X) 



- "\ f ''v.^' t = Pt\Born). (95) 

Thus Pt\0{as)] X) = Pi'^\Born) as long as one can neglect non-Born term like structures 
in the radiative as-corrections resulting either from the one-loop or the A-dependent hard 
gluon corrections. As it turns out the non-Born term like corrections are in general 
small but can amount to several percent. The above reasoning breaks down when either 
the numerator or the denominator in Eq. (p5!l approaches zero which can happen for very 
small values of A. As has been argued before such small cut values are not acceptable from 
the physics point of view. 

In Fig. Owe show a plot of P^ as a function of A/Amax again for the three cm. energies 
a/s = 400, 500 and 1000 GeV where P^ is the longitudinal polarization of the top quark 
P^ = cr^/o". Note that again one has to separately integrate the numerator and denomi- 
nator over the gluon energy when calculating P^, i.e. P^{X) = a^{X)/a{X). As in Fig. [8] 
the radiative corrections and the dependence on A can be seen to be quite small. The 
longitudinal polarization P^ is largest for y/s = 1000 GeV. 

In order to highlight the size of the radiative corrections to P^ we define a fractional 
deviation of P^ from its Born term value for different cut-off values by writing 

.p,, ^ P'(X) - P'iBorn) 

where P^{X) is the value of P^ for the upper cut-off parameter A, i.e. in our above termi- 
nology P^{X) refers to the value of the observable in the soft region. Fig. [10] shows that 
close to A = the fractional deviations S{P^) tend to infinity because the denominator in 
P^(A) = cr^(A)/(j(A) go to zero, as mentioned before. Away from A ~ the dependence 
of S{P^) on the gluon cut A is not very pronounced except for the highest energy value 
^/s = 1000 GeV. The fractional deviation is largest for ^/s = 400 GeV. 



•^The forward-backward asymmetry Aps defined in Eq. (|94p is such a polarization-type observable with 

(m) 
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Figure 8: Dependence of the forward-backward asymmetry Aps on the (upper) cut-off 
A/Ajnax in the soft region (full line). Curves are shown for the three center-of-mass energies 
= 400, 500, and 1000 GeV. Also shown are the respective cut-off independent LO Born 
term contributions (horizontal dashed lines). 
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Figure 9: Dependence of the longitudinal polarization on the (upper) cut-off A/Amax in 
the soft region (full line). Curves are shown for the three center-of-mass energies ^/s = 400, 
500 and 1000 GeV. Also shown are the respective cut-off independent LO Born term 
contributions (horizontal dashed lines). 
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Figure 10: Dependence of the fractional deviation of the longitudinal polarization S{P^) on 
the (upper) cut-off A/Amax in the soft region. Curves are shown for the three center-of-mass 
energies ^/s = 400 (dotted line), 500 (dashed line) and 1000 GeV (full line). 

Also of interest are the values of a rate function in the hard gluon region. To this end 
we define a lower scaled gluon energy cut-off Aiower and integrate from Aiowcr to the upper 
limit Amax = (1 — 0/2- As before this is effectively done by subtraction, i.e. a^"^^{hard) = 
a^"^\Xaia.x) — a^'^\X) since we have not separately listed analytical formulas for the hard 
gluon rates. We then define a forward-backward asymmetry ApB{hard) and a longitudinal 
polarization P^{hard) in the hard region by writing 



Apsihard) 

and 



a{forward) — a{hackward) 
a {forward) + a {backward) 



(97) 

hard 



PHhard) = — (9J 

^ ^ - hard ^ 



a 

In Fig. [11] we show a plot of AFB{hard) as a function of A/Amax again for the three 
cm. energies i/i = 400, 500 and 1000 GeV. As the lower cut-off tends to zero AFB{hard) 
reaches values very close to those of AFB{soft) in Fig. [8] showing that the non-Born term 
structures in the a^-radiative corrections are not very significant. Only for larger cut-off 
values does one find significant deviation from the Born term values. For example, for 
A/Amax = 0.6 and y/s = 1000 GeV one has a 30% deviation from the Born term value. 

Fig. [T2] shows the same plot for the longitudinal polarization P^. Similar remarks apply 
as in the discussion of Aps^hard) except that the dependence on the lower cut-off is not 
as pronounced as in Fig. [TTl Marked deviations from the Born term values only set in at 
larger values of A. 
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Figure 11: Dependence of the forward-backward asymmetry Aps on A/Amax in the hard 
region where A denotes a lower cut-off. Curves are shown for the three center-of-mass 
energies ^/s = 400 (dotted hne), 500 (dashed hne), and 1000 GeV (fuU hne). The straight 
hnes indicate the Born term level results. 
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Figure 12: Dependence of the longitudinal polarization P'' on A/Amax in the hard region 
where A denotes a lower cut-off. Curves are shown for the three center-of-mass energies 
= 400 (dotted line), 500 (dashed line), and 1000 GeV (full line). The straight lines 
indicate the Born term level results. 
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7.2 NLO contributions to vanishing LO observables or structure 
functions 

It was pointed out already in Ref. [8] that the longitudinal polarization of the top quark 
produced from a longitudinally polarized gauge boson (7 and/or Z) denoted by P[ vanishes 
at the Born term level. P[ vanishes at the Born term level and also for the one-loop 
contribution due to the two facts that there are no second-class currents in the SM and 
that one is dealing with a two-body final state in these two cases. Technically this comes 
about since the contractions of the first class axial currents uj^'j^v and uq^-^^v with the 
longitudinal projector 63 (see Eq. ( 125!) ) vanish in the two-body case. In the Standard 
Model a nonvanishing value of the polarization is generated only at NLO (or higher 
orders) from real gluon bremsstrahlung. This NLO effect is quite small as can be seen 
from Fig. 2a in Ref. [H] which shows that rises from zero at threshold to —0.21% at 
^/s = 1000 GeV. 

A larger absolute value of P[ is obtained in the hard gluon region since P[ is an 0{as) 
effect. To this end we define the ratio 

Pf(hard) = ^ (99) 

^ a hard ^ ' 

where the hard gluon region is defined as in the beginning of this section. In Fig. [13] we 
show a plot of P[{hard) as a function of the scaled gluon energy cut-off where the cut-off 
parameter A now refers to a lower cut-off. It goes without saying that P[{hard) = in 
the soft-gluon or eikonal approximation since then cr£(as) oc a\^{Born) = in the soft- 
gluon or eikonal approximation. Fig. [13] shows that P[{hard) can become as large as —4% 
for y/s = 1000 GeV and A/Amax = 0.8. P[{hard) increases when the energy increases. 
P[{hard) goes to zero as A — since in this limit af^ is finite whereas a diverges. 

We mention that a nonvanishing contribution to P[ can also be obtained by adding an 
anomalous axial current to the usual SM first class top quark current structure. This will 
be discussed later on. 

There are two classes of relations among the structure functions Hl^"^^ at the two-body 
level. The first class of relations depends solely on the fact that one is dealing with a 
two-body final state at the Born term and one-loop level. There are four relations of this 
kind 

real part: Hi; = H^^ Hi = HI^ H^ = Hfj^ (100) 
imaginary part: Hp = Hff (101) 

The second class of relations depends on the two-body dynamics and on the fact that one 
has only first class currents in the SM. There are six relations of this kind. These are 

real part: Hi = Hi iff = H'/ = Hf Hf = (102) 
imaginary part: Hf = Hf H]^ = Hj^ (103) 
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Figure 13: Dependence of the longitudinal polarization from a longitudinal polarized gauge 
boson P[ on A/Amax in the hard gluon region where A denotes a lower cut-off. Curves are 
shown for the three center-of-mass energies y/s = 400 (dotted line), 500 (dashed line), and 
1000 GeV (full line). 

One can explicitly check with the Born term and one- loop expressions listed in Sec. 3 that 
these relations are in fact satisfied. 

Note that the class 1 relation Hp = H^f, and the class 2 relations = Hf^ and 
Hj^ = Hj^ will not be affected by the 0{as) tree graph contributions since they result 
from the imaginary parts of the (two-body) one-loop contributions. As mentioned before, 
the relation if^ = H"^ interestingly also holds at the 0{as) tree graph level. In the 
following we shall numerically investigate how the remaining relations in fllOOl) and fll02p 
are affected by the 0{as) tree graph contributions. It goes without saying that the relevant 
remaining relations in fllOOp and (11021) still hold at NLO if one uses the soft-gluon or eikonal 
approximations rather than the exact form of the radiative corrections. 

We start our numerical discussion with the first class of relations in Eq. fllOOl) . In order 
to obtain a quantitative handle on how the tree graph contributions affect the first class 
relations HIj = H]^ and Hp = H^f in Eq. fllOOl) we consider differences of the relevant 
structure functions and (arbitrarily) normalize them to Hlj{Born). In Fig. [THwe show a 
plot of the ratios {Hlj—Hp)/Hlf{Born) and {Hp—H^)/Hjj{Born) as functions of the upper 
cut-off ("soft region") in terms of the scaled gluon energy cut A/Amax for y/s = 500 GeV. 
The violation of the class 1 relations slowly rises from zero at the soft-gluon point and 
reaches values of 0.27 and —0.02%, respectively, for the two above ratios at Amax where one 
integrates over the full gluon phase-space. In Fig. [TS] we consider the hard region where 
A/Amax now refers to a lower cut-off in the gluon energy. Now A = corresponds to a 
full phase-space integration and one therefore recovers the A/Amax = 1 limiting values of 
Fig. [m remembering that there are no loop contributions to the above four quantities. The 
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relevant ratios go to zero for A = Amax in Fig. [15] since phase-space goes to zero. 

The influence of the tree graph contributions on the second class of relations Eq. (I102p is 
tested in a similar manner. We consider again differences of the relevant structure functions 
(or structure functions themselves) normalized to Hlj{Born). In Figs. fT6] and [T71 we show 
plots of the ratios {Hi - Hi) / Hl^Born), Hf / Hlj{Born), {Hf - Hf)/Hlj{Born), and 
[Hf^ — H^)/Hlj{Born) for upper and lower cut-off values of the gluon energy, respectively. 
In Fig. [16] ( "soft region" ) the violations rise from zero at the soft-gluon point to the values 
0.29, 0.15, 0.08, and 0.03% for A = Amax where one integrates over the full gluon phase- 
space. Fig. [T71 shows the same four ratios in the hard gluon region. As before the right-most 
values in Fig. [16] agree with their left-most pendants in Fig. [17] The violations of the class 
1 and class 2 relations due to hard gluon radiation can be seen to be generally quite small. 

The effect of the radiative corrections to the class 2 relations fll02p can be mimicked by 
adding an anomalous axial current to the SM currents. The anomalous axial current to be 
added reads (see e.g. pTl [TB] ) 

- ia^^qy 

f {anomalous) = Qaipt— — -Ibipt- (104) 

2mt 

In general Qa can be complex, Qa = Rega + ilmga- Note that the current in Eq. fll04p 
is a so-called second-class current with J^*" = 1^^ quantum numbers. In particular, the 
contraction of the anomalous current with the longitudinal projector 63 (see Eq. (j25l) ) no 
longer vanishes, i.e. one now has e'i,uafj_uq"v 7^ 0, and therefore Hjf 7^ 0. It should be 
clear that the addition of the anomalous axial current does not affect the class 1 two-body 
relations in Eq. fllOOp but, in general, violates the class 2 relations. We assume that the 
coupling strength ga is small and we therefore only consider the interference contribution 
of Eq. fll04p with the SM {tt) current, i.e. terms that are linear in ga- 

The interference contribution of the anomalous axial-vector current can be calculated 
using the projection formulas written down in Sec. 3. One finds 

Hi- Hi = 0{gl) 

Hi' = -2N,q\Rega 

H^-Hf = N.^qWRega 
y 2m 

Hf-Hf = N^^q\Rega. (105) 
V 2m 

It is noteworthy that only the real part of ga contributes to the relations fll05p . In 
order to obtain a quantitative handle on the coupling parameter ga we determine the 
values of the anomalous parameter ga that would reproduce the fully integrated quanti- 
ties Hl'/Hlj{Born), {Hf - Hf) / Hlj{Born) and {Hf - Hf)/Hlj{Born), i.e. the values 
that these quantities take at the right-hand side of Fig. [16] at A/Amax = 1- One finds 
ga = -0.0032, 0.0023 and 0.0007 for i/f , H^ - Hf and Hf - Hf^, respectively. Values 
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Figure 14: Goodness of the class 1 relations against radiative corrections using an upper 
gluon energy cut. Dependence of the ratios (iJ^ — Hp)/ Hlj{Born) (sohd hue), and {Hp — 
H^f)/ Hlj{Born) (dashed line) on X/Xmax where A denotes an upper cut-off. Curves are 
shown for the center-of-mass energy y/s = 500 GeV in the soft region. 
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Figure 15: The same as in Fig. [T4]for the hard region. 
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Figure 16: Goodness of the class 2 relations against radiative corrections using an up- 
per gluon energy cut. Dependence of the ratios {H}^ — H\) / Hij{Born) (solid line) 
Hf/Hlj{Born) (dashed line) {HY - Hf)/Hlj{Born) (dash-dotted line), and {Hf - 
H^) / Hlj{Born) (dotted line) on A/Amax where A denotes an upper cut-off. Curves are 
shown for the center-of-mass energy y/s = 500 GeV. 
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Figure 17: The same as in Fig. [16] for the hard region. 
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substantially larger than these combinations of structure functions would signal contri- 
butions from a second-class current with coupling strength exceeding the above values of 

9 a- 



8 Summary and outlook 

We have presented analytical results for the 0{as) radiative corrections to polarized top 
quark pair production in e^e~ annihilation with a specific gluon energy cut. When the 
gluon energy cut is taken to its maximal value we recover previously known results [3, [TT] . 
The size of the radiative corrections to polarization-type observables involving the top 
quark is generally quite small in the soft-gluon region but can become substantial in the 
hard gluon region. This in turn implies that the dependence of the polarization-type 
observables on the gluon energy cut is generally quite small in the soft-gluon region but 
can become large in the hard gluon region. We have calculated the contributions of a CP- 
odd non-SM coupling to some linear combinations of structure functions that vanish in 
the two-body SM case. These were compared to SM contributions resulting from radiative 
corrections. 

We have not considered beam polarization effects in our analysis. However, in as much 
as we have calculated the complete set of single spin structure functions, beam polarization 
effects can be easily incorporated into our analysis as described e.g. in more detail in Ref. [8]. 

We have decomposed the top spin vector in the helicity basis, i.e. the z direction of our 
spin basis is determined by the momentum of the top quark. In addition to the helicity 
basis the authors of Refs. [161 [IH] have also considered a beamline and an off-diagonal 
basis. A discussion of how these bases are related to the helicity basis in the context of 
the NLO corrections can be found in Ref. [TT] . 

All the results in this paper refer to the polarization of the top quark. In order to obtain 
the SM and anomalous coupling predictions for the polarization of the antitop quark let 
us first set up an orthonormal spin basis for the antitop quark by replacing the momenta 
in Eq. (jlj) by their charge conjugate partners, i.e. pi — p2 and p^- — > Pe+- The three 
orthonormal basis vectors (cr, e/v, e^) are now given by 



In the polar angle distribution Eq. ([8]) the polar angle now refers to Oie- and not to 6 = 6te- 
as in the top quark case discussed in the main part of this paper. Since the lepton pair is 
back-to-back, one has Off,- = vr — 6te+, i-e. the two terms in Eq. proportional to cos 6 
change sign if written in terms of cos 6te+ - 

Let us list the SM Born term and the anomalous contributions in the antitop quark 
case given by Eq. (11041) together with the relevant contributions in the top quark case. 
One finds 




Pe+ X P2 
\Pe+ X P2 




(106) 
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H^' = 0, Hlf = ±ANcq'v, 

^9a\ rT2 AT .2/1 ..2\ at „.2 2 IS'al 



= Ncq\l - v') + Ncq'v'^, Hi = Ncq\l - v') - Ncqh 

Hi' = -2Ncq'vImga, Hi' = -2Ncq'vReg,, 

Hi' = 2Ncq\l + v^), Hi' = 2Ncq\l - v"), 

Hi = 0, Hl = ±^N,q\, 

Hi = —^'vlmga, Hj = — =w(±^ + Rega) 



2i ^ ^ 

H7 = ± v'Reg^), Hf = ^{iT v'Reg^), 



2^ ^ ^ 

TjlN t^Cq 2t tt2N ^^Cq 2t 

Hj = ±^=v Imga, Hj = T^=^v Imga, 



Hr = ^vi±^-Reg^), Hf = ^vlmg^, (107) 



where the upper and lower signs refer to the top quark and antitop quark cases, respectively. 
As concerns the SM Born term contributions one finds 



at (cos 9te-) = cif (cos 9te+ ) 

Pt'^ {cos 9 te-) = -PI'^ {cos 9te+) 

P^{cos9t,-) = P^{cos9te+). (lOJ 



In the three-body case one has to simultaneously exchange {y ^ z) in the SM part of 
Eqs. fll07p and (11081) . For example, one has hI''^' {top; y , z) = —HI''^' {antitop; z,y). If one 
performs an integration symmetric in y and z as done in this paper the SM part of the 
relations (11071) and (llOSp also hold for the integrated three-body results. 

The linear contributions of the anomalous coupling to the polarization vector behave 
in the opposite way to those in Eq. (11081) . i.e. 

P''^ {anomalous; cos 9 te') = P^^ {anomalous; cos 9 te+) 
Pj {anomalous; cos 9 te-) = —Pf {anomalous; cos 9 (109) 

It is clear that one can obtain an additional handle on the anomalous contributions by 
taking sums and differences of the top quark and antitop quark polarizations. For example, 
{Pf'^ {cos9te-) + -Pf '^(cos 6'fe+)) and {P^ {cos9te-) — Pf{cos9te+)) are contributed to only 
by the anomalous contributions. 

In this paper we have not discussed how the spin of the top quark can be analyzed. The 
top quark decays weakly and is therefore self-analysing. If one assumes SM interactions 
in the cascade decay t — > bW^{^ l^i^i, qq) the polarization of the top quark can be 
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reconstructed by measuring spin-momentum correlations either in the top quark rest system 
(see e.g. Refs. [20l[2Tl[22l[23]) or in the W rest system as e.g. discussed in Refs. [T5l [2it [25] . 
We mention that there exists a large body of literature of how non-SM interactions in the 
production (see e.g. Ref. [26j) (such as the anomalous coupling Eq. fll04p ). and/or in the 
decay affect such spin-momentum correlations (see e.g. Ref. [27] and references therein). 

Gluons can be emitted from the original production process e"*"e~ —>■ ti{G) as well as 
from the follow-up decay process t ^ b + W~^{G) and t + W~{G) where we take the 
W's to decay leptonically. Interference effects between the two processes are expected to 
be quite small since they are suppressed by a factor of ~ Tt/rrit ~ 1%. In order to identify 
the gluons of the original production process (which are the subject of this paper) one has 
to demand that the gluon's four-momentum satisfies q = Pt + Pt + Pg- Gluons that satisfy 
Pt = Pb + Pw + Pg or Pt = Pb+ Pw + Pg clearly originate from the follow-up processes and 
can thus be vetoed. How effectively gluons not originating from the original production 
process can be removed from the data sample has to be carefully studied in detailed Monte 
Carlo simulation runs. 

With the appropriate modifications our results can also be apphed to the (bb) case. 
While the Imxz contributions resulting from the imaginary part of the Breit-Wigner line 
shape are negligibly small in the [tt) case (since (tt) threshold is far away from the Z pole) 
the Imxz contribution is more pronounced in the (bb) case in particular in the vicinity of 
the Z pole. However, close to the Z pole the transverse and normal polarization of the 
bottom quark are severely suppressed due to the overall helicity suppression factor 2m/ ^i. 
In this sense the phenomenology of the top quark spin above (ti) threshold is richer than 
that of the bottom quark in the high energy realm. 
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A SM values of the electroweak coupling coefficients 

The electroweak coupling matrix elements gij{q^) are given by 

911 = Q}-2QfVeVfRexz + {vl + al){vj + a})\xz\'', 

912 = Q}-2QfVeVfRexz + {vl + al){vj-a})\xz\'^, 

913 = -2QfVeafImXz, (Al) 
9u = 2QfVeafReXz -2{vl + al)vfaf\xz\^: 



921 = Q'f-2QfVeVfRexz + {vl-al){vj + a'j)\xz 

922 = Q'f-2QfVeVfRexz + (v^-al)(vj-aj)\xz 

923 = -2QfVeafImXz, 

924 = 2QfVeafReXz-2{vl-al)vfaf\xzf, 



931 = -2QfaeVfImXz, 

932 = -2QfaeVfImXz, 

933 = 2QfaeafReXz, 

934 = 2QfaeafImXz: 



(Al) 



5'42 = 2QfaeVfReXz -2veae{vj - a'j)\xz\'^, 

943 = 2QfaeafImXz, 

9u = -2(5/aea/i?eXz + 4veaeV/a/|xz|^ 

where Xz(?^) = gM^q^ /{q^ — M| + zM^F^), with Mz and the mass and width of the 
Z° and g = G'F(8\/27ra)~^ ?a 4.49- lO^^GeV"^. Q/ are the charges of the final state quarks 
to which the electroweak currents directly couple; fg and ag, Vf and a/ are the electroweak 
vector and axial-vector coupling constants. For example, in the Weinberg-Salam model, 
one has Vg — —1 + 4sin^ 9w, = ~1 for leptons, Vf — 1 — ^ sin^ 9w, = 1 for up-type 
quarks {Qf — |), and Vf — —1 + ^ sin^ 9^, o/ = —1 for down-type quarks {Qf — — |)- The 
left- and right-handed coupling constants are then given by 9l — v + a and 9r = v — a, 
respectively. In the purely electromagnetic case one has gn = gi2 = g2i = g22 = Q} and 
all other g^'r = 0. The terms linear in Rexz and ImXz come from 7 — interference, 
whereas the terms proportional to |xzp originate from Z exchange. 

Contributions coming from the imaginary part of the Breit-Wigner resonance shape 
arc of order 0{ImXz{q^) I ReXz^q^)) and can thus safely be neglected for top quark pair 
production. For example, in the threshold region of top quark pair production ImXz/ReXz 
is approximately 0.1% and decreases further with a. 1/q^ power fall-off behaviour. 
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B Decay rate terms ti 

It is convenient to define the mass dependent variables a := 2 + v^j b := 2 — y/^ and 
w := ^(1 — \/^) /(I + \/^). The rate functions ti, . . . ,ti2 appearing in the main text are 
then given by 

" - Ks^)- '--KtTT!) - (I) (^^) 

(2 CI 

t4 := Li2(w;) - Li2(-w) + Li2(-^t«) - Li2(--u') (B3) 

* 2 l^4(l + Vf)/ Vl-t>; 'Va(l+jf)y 'l^a(l-ti;)y 

^L, (l±i^) - Li, (i^) + L. - L, (^) (B4) 

tQ := ln2(l + w) + ln^(l - w) + In ln(l - w^) + 

+Li,(^)+Li,(^)-2Li,(^) (B5) 
21, I In ri±^^ - Li. (^] + L., f-^^ + 



/ 2aw \ ^ / 2aw \ ^ / aw \ „^ / aw 
+Li2 - Li2 - 2Li2 + 2Li2 



b + awj \ b — awj \b + aw/ \ b — aw 

- m (lU (l±^) + Li. (J^) - Li, (- A) - (BT) 



+3rLi,f-^)-Li,f^)) (B8) 
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C Decay rate terms £i, to±, ti±, and t 



w 



The logarithmic rate terms ii arc given by 

1 + till \ 



^2 



14+ 



4+ 
4+ 
4- 
4+ 
4+ 



In 



'wj — wl' 



= In 



wt 



In 



+ In 



awi 
aw2 



I - W2 



— In 
-In 



1 - 2A + ve; 



In 



(-) 

\WiJ 

AC 

■— + — + 2v 

yi y2 



1 + V 



Wo - Wi, 



1 -t; 







In(^) +21n. 2 



Wo - W2J 
- 1 



+ 



2v 



2-21„(H^)+ln 



(1 + wi)(6- awi)_\ _^ aw2)(l - wa)^ 



= In 



+ (2^;-(2-e)ln 

1-1^2 



1 



V 



Wq 



In 



21n(^^) - In 
1 — V 



1 - V, 
1 + Wi 

1 + Wo 

1 + Wi 

6 — awi 



Wq — Wf 



4- = 21n 



wg - wi 



Wn 



wr 



1 

1 + v 



Wi 



- 1 



-In I 

F aw2 



' b + aw2 ' 

1 - W2 , 



In 



= In 



'w^-wl 



wj 



Wo - W; 



i9+ = In 



^Wq — Wi^ 

'l-w|V 



In 
In 



1 



Wi 



i^- = In 



■W^2 / 
W2 - Wa 
Wi - Wx 
Wq - W2 
Wq - Wi 

1 - W2 

1 — Wi 



-In 

- In 



W2 + Wa \ 

wi + Wa/ 

/ Wq + W2 \ 

Vwq + Wi/ 



(CI) 
(C2) 
(C3) 
(C4) 



1 + wJ 



+ 

(C5) 
(C6) 
(C7) 
(C8) 
(C9) 
(CIO) 
(Cll) 



while for the additional phase-space contribution we have to use 



•''2 



lnfl±^Ul„f^±^') 
\1 — W2/ \b — aw2j 

il_ = ln(l-w^)+ln(6'-a'w^), 

il_ = ln6 = ln(2-ye), 

l%_ = ln(l -w^) -ln(6' -a'w^), 



£1+ = In 



Wq + W2 
Wo - W2 



W2x 

In 



Wn 



Wq - W| 



(C12) 
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For the double and dilogarithmic decay rate terms we obtain 







- t'^{w,) - 1\ 








^Oit — 










-CM)- 


tl± = 













(C14) 



(1 + V^)V^\ fWo±W2 



iba 



+ + (C15) 



while for the additional phase-space contribution we take 

C = t"(^2) - e(0), = C±(^2) - C±(0), = ^^±(^2) - iliO) (C16) 

where 

t^iw) = Li2H-Li2(-^) + Li,(^)-Li2(^), 



= 21n(wo) M'f^) + Li2(i^') - Li2(-u.') - Lis (^;^) - Li2 



C» = 21n(6)lnH + Li2H-Li2(-«;)-Li2('^') -Li/ """^ 



v) 



ti{w) = ln^H + Li2H + Li2(-w) + Li2(^)+Li2(^), (C17) 



ttH = -2 In (^^) ln(wo - w) + 



+Li2 r^^^ - Li2 C^^) + Li2 f - Li2 " 



\wo + 1 J \ Wo — 1 J \ awo + b J \ awo — b 

'1 + 



= -2 In ln(wo + w) + 



+Li2 f - Li2 f ^^^1 + Li2 f + - Li2 + " 



\wo + I J \ Wo — 1 J \ awo + b J \ awo — b 

47 



to_{w) = ^ In (j^-^ ln(^o ln^(^o -w) + 



2wo J \ Wo — 1 J Vwo + 1 

to^w) = ~\n(^j^^\n{wo + w) + ^\n\wo + w) + 



-Li 



/ V Wo - 1 / V Wo + 1 



2w, 

t^t{w) = -InCHwo-w) + 

+Li, f "^^"-^^ "l + Li, f "(""-^^ l - Li, C-^^) - Li, f^^) , 
y awo — / V "'"^0 + Vwo — 1/ Vwo + l/ 

-Li, ( + _ Li, f + + Li, f + Li, (^^^) , 
V awo — J \ awo + o I \Wo — lJ Kwq + IJ 



+Li, ( '"'-'" ) + Li, ( - Li, (^^) - Li, f ^^^i^' 

to+{w) = W(wo + w) - ln(l - Wq) ln(wo + w) + 

+Li, (1^5L^) + Li, (!^) + Li, (^^5L±^) + Li, (!^) , (C18) 

\Wo + W J \Wo + W J \Wo — lJ \Wo + l / 

tl''_{w) = ln'(l-w) + lnQ)ln(l-w) + 

i?^(w) = ln2(l + w) + In ln(l + w) + 

+Li,f4^ULi,f^(i±^ULi,^^ + '" 



a(l + V 4 y ' V 2 



^ -,n(i±-)K:-.)-U,(i^).Li,(i^)-L.,(^ 
. ,„(i±-),,,.)_Li,(i±ii),Li,(i±^).Li,( 
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a{l-w) J \ 4 J " V 2 



t-W - K2a).n(l + ») + Li,(^)-L.(^)+L,(i±i^;) 

^^_(^/;) = ^ln2(l --«;) +ln21n(l -w) -ln(l -'u;^)ln(l --«;) + 
^. , l-w\ ^. fl-w\ ^. fl-w 



tl^{w) = -ln2(l +W) - ln21n(l + w) + 
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